Attosecond electron thermalization by laser-driven electron recoUision in atoms 
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Nonsequential multiple ionization of atoms in intense laser fields is initiated by a recoUision 
between an electron, freed by tunneling, and its parent ion. Following recoUision, the initial electron 
shares its energy with several bound electrons. We use a classical model based on rapid electron 
thermalization to interpret recent experiments. For neon, good agreement with the available data 
is obtained with an upper bound of 460 attoseconds for the thermalization time. 
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PACS numbers: 

Atoms exposed to intense laser fields ionize. The freed 
electron and its ionic partner are accelerated by the laser 
field away from each other. When the field changes sign, 
it may drive the electron into a recoUision with the ion. 
This simple mechanism, which for high-intensity low- 
frequency fields is largely classical, governs many laser- 
atom processes such as high-order harmonic generation 
(HHG), high-order above-threshold ionization (HATI), 
and nonsequential double ionization (NSDI) and explains 
the gross features of the spectra observed _lj . 

The period of the commonly applied titanium-sapphire 
(TiiSa) laser is about 2.7 fs. The recoUision physics un- 
fold on the time scale of a small fraction of the laser pe- 
riod. Therefore, the analysis of laser-induced recoUision 
phenomena provides access to the inner-atomic dynamics 
on the attosecond time scale and, indeed, the focus of re- 
cent investigations has moved into this temporal domain. 
For example, it has brought molecular imaging with sub- 
angstrom spatial and subfemtosecond temporal resolu- 
tion within reach ^2, .3i]. The advent of phase-stabilized 
infrared few-cycle pulses and of uv pulses of attosecond 
duration allows even more control but neither are 
necessary for a study of the attosecond dynamics. 

In this Letter, we analyze recent experiments on nonse- 
quential multiple ionization (NSMI) of neon [5,6] in which 
the momentum distributions of Ne^+ [N — 3, 4) were 
measured at two intensities. We use the fact that the 
time-dependent laser field provides a clock — the field 
accelerates the ion to a final velocity that depends on 
the times at which the N electrons ionized. This "streak 
camera" |2| therefore measures the range of times of ion- 
ization. Comparing the measured momentum distribu- 
tions to those predicted by a classical model, we infer that 
the recoUiding electron thermalizes with the iV — 1 bound 
electrons in less than 500 attoseconds. To our knowledge, 
no other method allows such a low upper bound for ther- 
malization times within atoms to be measured. 

Nonsequential double and multiple ionization is de- 
fined by the fact that it is not sequential, that is, it is 
not the product of a sequence of uncorrelated single- 



ionization events. NSDI and NSMI require electron- 
electron correlation as a necessary precondition §]■ 
Even for the very simplest such process - NSDI of helium 
- a fully quantum-mechanical description from first prin- 
ciples, i.e. by solution of the time-dependent Schrodinger 
equation in six spatial dimensions, has not been accom- 
plished yet '9], and for the heavier atoms it is clearly 
out of the question. This leaves approximate quantum- 
mechanical approaches, such as density-functional meth- 
ods |lOl |. 5'-matrix methods that try to identify the most 
relevant terms of an appropriate perturbative expansion 
[lllIT^ . or classical-trajectory methods T^Q]. For mul- 
tiple (triple and higher) ionization, which occurs under 
the same conditions as NSDI if the laser intensity is high 
enough, any description from first principles appears to 
be utterly out of reach. It seems equally hopeless, for 
three electrons and more, to identify the relevant dia- 
grams in the microscopic S'-matrix approach. 

The model we here propose is in the spirit of the one 
of Ref. 0, but essentially classical. We investigate the 
scenario wherein NSMI is effected by one single recoUi- 
sion. We assume that the pertinent electron tunnels into 
the continuum with zero velocity at the ionization time t' 
according to the time-dependent v&teRit'), for which we 
adopt the standard quasi-static rate [l^- Thereafter, we 
turn to an entirely classical description: The laser field 
may drive the electron back to its parent ion at a later 
time t, which is a function t{t') of the ionization time 
and can be easily evaluated. We assume that the en- 
ergy i?i-et {t) of the returning electron be completely ther- 
malized among the ensemble of participating electrons, 
that is, the returning electron and the — 1 electrons to 
be freed. These A^ electrons then form an excited com- 
plex with the total energy (with respect to the continuum 



threshold) i?ret(i) — £'o^\ where E'^" ' > denotes the 
total ionization potential of the A^ — 1 (up to the recol- 
lision time t inactive) electrons. The distribution of en- 
ergy and momentum over the A^ electrons is assumed to 
be completely statistical and only governed by the avail- 
able phase space. At the time t -\- At, the A^ electrons 
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become free to move in the laser field, which is described 
by the vector potential A{t) such that A{t) = outside 
the pulse. 

The corresponding distribution of the final electron 
momenta p„ (n = 1, . . . , TV) is proportional to 

F(pi,p2,...,PAr) = J dt'R{t')S(^E^o''^ -E,,t{t) 

1 ^ \ 

+ -^[p„ + A(t + A<)p , (1) 
ji=i / 

where the integral extends over the ionization time t' . 
The S function expresses the fact that the total kinetic 
energy of the N participating electrons at the time t + At 
is fixed by the first-ionized electron at its recollision time 
t. This constitutes the one and only condition on the final 
momenta p„. The only free parameter of this model is 
the time delay At between the recollision time and the 
time when the electrons become free. It is the sum of a 
thermalization time Atth ~ the time it takes to establish 
the statistical ensemble - and a possible additional "dwell 
time" , until the electrons become free. By comparing the 
predictions of the model with the data, we will be able 
to infer a value of At, which in turn provides an upper 
limit for the thermalization time Atth- 

This model is an extension to NSMI of a classical 
model introduced for NSDI for At = [3 E3- Suf- 
ficiently high above threshold, it produced momentum 
distributions that were virtually indistinguishable from 
their quantum-mechanical counterparts. Statistical mod- 
els similar to the one above have been used in many areas 
of physics. For example, the statistical Rice-Ramsperger- 
Kassel-Marcus (RRKM) theory ^ describes thermal- 
ization of molecular vibrational degrees of freedom, and 
for high-energy collisions of elementary particles and 
heavy ions statistical models have been utilized to predict 
the momentum spectra of the reaction products [l9j . An 
excited complex as the doorway to NSMI was also con- 
sidered in Ref. t2Qi] . 

A convenient feature of the ansatz is that integra- 
tion over unobserved momentum components is easily 
carried out. To this end, we exponentialize the S func- 
tion in Eq. with the help of its Fourier representation 



S{x) = [ — exp(— iAx). 



(2) 



Infinite integrations over the momenta p„ can then be 
done by Gaussian quadrature. The remaining integration 
over the variable A is taken care of by the formula |2ll | 
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where a;"^ — x"9{x), with 9{x) the unit step function and 
+0. 



For comparison with the experiments 0, 0, we cal- 
culate the distribution of the momentum P of the ion. 
Provided the momentum of the absorbed laser photons 
can be neglected, momentum conservation implies P = 
-I]n=iPn' so that 

Fion(P)= j Wd^'VnS fp + ^p„^ F(pi,p2,...,pw) 
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dt'R{t') (A^;jv,ion)^~n4) 



with A^ATaon = E,,,{t) - E^'^ - - ^Mt + At)]2. 

In the experiments thus far, the ion momentum trans- 
verse to the direction of the laser polarization is entirely 
or partly integrated over. In the first case, the remaining 
distribution of the longitudinal ion momentum P|| is 



(27r) 



AF((3A- l)/2 



Eon{P\\) = j d^P^Eon(P) 

dt'Rit') (Ai?i„„||)J~^ 



(5) 



where now AE\, 



Eret{t) — E^ 



(AT) 



2-7v[^ll 



NA{t- 



At)]^. If just one transverse- momentum component 
{P^_,2^ say) is integrated while the other one (-Pi.i = P\_) 
is observed, the corresponding distribution is 



(27r) 



AF(3iV/2- 1) 



dt'R{t') {AEn^ 



ion||±j 



(6) 



^Ar,ion||_L = A£'ion|| " 2N^\ 



with AEjs 

In Fig. ^ we present calculations of the ion-momentum 
distributions for triple (upper panel) and quadruple 
(lower panel) NSMI of neon according to Eq. jSJ, for 
various values of At between and 0.2T. The parame- 
ters are for the experimental data of neon presented in 
Fig. 2 of Ref. |^. The figure displays the double- hump 
structure of the ion momentum in NSMI of neon. When 
the time delay At increases from 0, initially, the effect on 
the momentum distribution is small. Later, however, the 
center positions of the humps start moving towards zero 
momentum and the widths of the humps increase until 
the two humps begin to merge. This behavior can easily 
be understood from the recollision kinematics, which are 
illustrated in Fig. [21 The final electrons undergo maxi- 
mal acceleration if they are released near a zero crossing 
of the electric field, which in the figure occurs at t = T. 
For the example of triple ionization of neon, the earliest 

(3) 

recollision with E^ct > Eq takes place at t = 0.74T. 
Already with a delay of At > 0.26T, all electrons will be 
released after the zero crossing, which results in signifi- 
cantly lower ion momenta. 
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FIG. 1: (Color online) Distribution of the longitudinal ion mo- 
mentum for triple (upper panel) and quadruple (lower panel) 
nonsequential ionization of neon at 2 PWcm"^ calculated 
from Eq. lO for various delays At as indicated in the lower 
panel. Note that in the upper panel the curves for At — 
and At — O.IT almost completely overlap. 



0.41 T 




t'/T 

FIG. 2: The kinetic energy Erct{t) (lower panel) and recol- 
lision time t (upper panel) of electrons liberated at the tun- 
neling time t' . The classically allowed intervals Atrcc of the 
recoUision time, for which the corresponding electron return 
energy Erct{t) is greater than Eq^^ (and -Eq*'), are shown in 
the upper panel. The calculation is for neon at 2.0PWcm~^. 



In Fig. 121 we compare the results of the statistical 
model with the data of Ref. . The intensities are those 
given for the experiment. We display momentum dis- 
tributions calculated from Eq. jSJ for zero delay and for 
At = 0.17T. The latter value was chosen to yield optimal 
agreement for the entire set of data. We notice, in partic- 
ular, that with this nonzero delay the model reproduces 
the maxima of the experimental ion-momentum distribu- 
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FIG. 3; (Color online) Distribution of the longitudinal mo- 
mentum of nonsequential triple ionization of neon at 1.5 
PWcm^^ (upper panel), 2.0 PWcm~^ (middle panel), and 
of nonsequential quadruple ionization of Ne at 2.0 PWcm~^ 
(lower panel). The rugged (black) curves represent the data of 
Fig. 2 of Ref. 0. The outermost smooth (red) curve and the 
innermost smooth (green) curve are calculated from Eq. Q 
for At = and At = 0.17T. 
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FIG. 4: Ion-momentum distribution of nonsequential triple 
ionization of neon at 1.0 PWcm~^ for At = (upper panel) 
and At — 0.17T (lower panel), calculated from Eq. ®, to be 
compared with the data of Ref. Q . 
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tion. This removes a longstanding discr epancy between 
models of the type discussed in Refs. [l^, lla UM ^^'^ 
data. For triple ionization, the calculated momentum 
distributions are wider than those of the data, in partic- 
ular for the higher intensity, even though the data may 
include a contribution of the partly sequential channel 
Ne — > Ne"*" Ne'^"'", or the contribution of a recoUision- 
excitation channel, both of which are, of course, not part 
of the model. 

Figure 0] exhibits the results of our thermalization 
model for = 3 for the conditions of Ref . . The dis- 
tribution of two components of the ion momentum is pre- 
sented, the component parallel to the laser field and one 
transverse component, which therefore provides a more 
stringent test of the model. The third component is in- 
tegrated over in the data, which corresponds to Eq. ((HJ. 
The intensity given in the experiment is 1.5PWcm~^; 
we obtain good agreement with the data for the lower 
intensity of l.OPWcm-^ and At = O.ITT as before 
The nonzero delay At has little effect on the transverse 
width of the distribution, but it causes an elongation in 
the longitudinal direction and moves the centers to lower 
momenta, markedly improving the agreement with the 
data. 

Encouraged by the good agreement between the model 
and the data, we interpret the delay for which we ob- 
served optimal agreement as an upper bound of the 
thermalization time Atth, as argued above. This yields 
A^th < Atopt = O.ITT ~ 460 as. A lower hmit of the 
thermalization time should be given by the inverse of 
the plasma frequency for an electron density pc — N in 
atomic units. This produces Atth < y'Tr/N, which is of 
the order of the atomic unit of time. 

All of the above discussion has been for neon. NSDI of 
argon appears to be governed by a recoUision-excitation 
scenario |23|. We note in passing that for (quadruple) 
NSMI of argon we get good agreement of our statistical 
thermalization model with the data Q for the experi- 
mental intensity and At — 0.35T, twice as long as for 
neon. Details will be given elsewhere. 

To test the model further, and possibly to set a tighter 
upper limit on the thermalization time, it is necessary to 
restrict the time range A^cm of recoUision (Fig. I^J. If 
thermalization is more rapid than 460 as, then the width 
of the ion-momentum distribution will decrease as we re- 
strict Aircc- There are a number of ways to minimize 
A^rec- Within limits, all that is needed is to lower the 
light intensity or to increase the laser frequency. How- 
ever, the best way to minimize Atj-cc and to control the 
time of recollision is to use a second harmonic field, po- 
larized perpendicular to the fundamental. The combined 
requirement that the electron and ion recoUide in both 
directions allows the time of recollision to be precisely 
determined and controlled via the relative phase of the 
two beams. 

The limit on the thermalization time that we deter- 



mine (as well as the much tighter bounds that seem fea- 
sible in future experiments) should apply to stationary 
electron-atom scattering in general. From a collision- 
physics perspective, as a result of streaking, the time- 
dependent laser field reveals information that is hard to 
obtain by other means. The streaking principle should 
also be applicable for studying nuclear dynamics. As in 
NSMI, nuclear processes also can be initiated by laser- 
controlled re-collision ji^l (of course, at much higher in- 
tensity) . Any nuclear decay process that results in a mass 
or charge change of the fragments will be streaked by the 
laser field just as electrons and ions are streaked in our 
case. 

To summarize, by comparison of experimental data 
for triple and quadruple nonsequential ionization of neon 
with a simple statistical recollision model where the re- 
turning electron thermalizes with a subset of the bound 
electrons, we have been able to conclude that (i) for 
neon such a model appears to contain the most relevant 
physics, and (ii) the time for this thermalization to occur 
is extremely fast, well below one femtosecond. 

We gratefully acknowledge stimulating discussions 
with G.G. Paulus, H. Rottke, and W. Sandner. 



[1] P.B. Corkum, Phys. Rev. Lett. 71, 1994 (1993). 

[2] H. Niikura, F. Legare, R. Hasbani, A.D. Bandrauk, 
M.Yu. Ivanov, D.M. Villeneuve, and P.B. Corkum, Na- 
ture (London) 417, 917 (2002). 

[3] S.N. Yurchenko, S. Patchkovskii, I.V. Litvinyuk, P.B. 
Corkum, and G.L. Yudin, Phys. Rev. Lett. 93, 223003 
(2004). 

[4] A. Baltuska, Th. Udem, M. Uiberacker, M. Hentschel, E. 

Goulielmakis, Ch. Gohle, R. Holzwarth, V.S. Yakovlev, 

A. Scrinzi, T.W. Hansch, and F. Krausz, Nature (Lon- 
don) 421, 611 (2003). 
[5] R. Moshammer, B. Feuerstein, W. Schmitt, A. Dorn, 

CD. Schroter, J. Ullrich, H. Rottke, C. Trump, M. 

Wittmann, G. Korn, K. Hoffmann, and W. Sandner, 

Phys. Rev. Lett. 84, 447 (2000). 
[6] A. Rudenko, K. Zrost, B. Feuerstein, V.L.B. de Jesus, 

CD. Schroter, R. Moshammer, and J. Ullrich, Phys. Rev. 

Lett. 93, 253001 (2004). 
[7] E. Constant, V.D. Taranukhin, A. Stolow, and P.B. 

Corkum, Phys. Rev. A 56, 3870 (1997). 
[8] Th. Weber, M. Weckenbrock, A. Staudte, L. Spielberger, 

O. Jagutzki, V. Mergel, F. Afaneh, G. Urbasch, M. 

VoUmer, H. Giessen, and R. Dorner, Phys. Rev. Lett. 

84, 444 (2000). 
[9] K.T. Taylor, J.S. Parker, K.J. Meharg, and D. Dundas, 

Eur. Phys. J. D 26, 67 (2003). 
[10] D. Bauer and F. Ceccherini, Opt. Express 8, 377 (2001). 
[11] A. Becker and F.H.M. Faisal, Phys. Rev. Lett. 84, 3546 

(2000); ibid. 89, 193003 (2002). 
[12] R. Kopold, W. Becker, H. Rottke, and W. Sandner, Phys. 

Rev. Lett. 85, 3781 (2000). 
[13] L.B. Fu, J. Liu, S.G. Chen, Phys. Rev. A 65, 

021406(R)(2002); J. Chen and CH. Nam, Phys. Rev. 



5 



A 66, 053415 (2002). 

[14] P.J. Ho, R. Panfili, S.L. Haan, and J.H. Eberly, Phys. 
Rev. Lett. 94, 093002 (2005). 

[15] L.D. Landau and E.M. Lifshitz, Quantum Mechan- 
ics (Nonrelativistic Theory) (Pergamon Press, Oxford, 
1977). 

[16] C. Figueira de Morisson Faria, X. Liu, W. Becker, and 

H. Schomerus, Phys. Rev. A 69, 012402(R)(2004). 
[17] C. Figueira de Morisson Faria, H. Schomerus, X. Liu, and 

W. Becker, Phys. Rev. A 69, 043405 (2004). 
[18] See, e.g., W. Forst, Theory of Unimolecular Reactions 

(Academic, New York, 1973); P.J. Robinson and K.A. 

Holbrook, Unimolecular Reactions (Wiley-Interscience, 

New York, 1972). 



[19] R. Hagedorn, Nuovo Cimento Suppl. 3, 147 (1965). 

[20] K. Sacha and B. Eckhardt, Phys. Rev. A 64, 053401 
(2001); J. Phys. B 36, 3923 (2003). 

[21] I.S. Gradsteyn and I.M. Ryzhik, Table of Integrals, Se- 
ries, and Products (Academic Press, New York, 1980). 

[22] A delay of the same order has also been observed in 
classical-trajectory simulations of NSDI of helium ,.14^]. 

[23] V.L.B. de Jesus, B. Feuerstein, K. Zrost, D. Fischer, A. 
Rudenko, F. Afaneh, CD. Schroter, R. Moshammer, and 
J. Ullrich, J. Phys. B 37, L161 (2004). 

[24] N. Milosevic, P.B. Corkum, and Th. Brabec, Phys. Rev. 
Lett. 92, 013002 (2004); S. Chelkowski, A.D. Bandrauk, 
and P.B. Corkum, Phys. Rev. Lett. 93, 083602 (2004). 



